General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



NASA 

Technical Memorandum 80251 

1HASA-TM-802S1) CRUSTAL EVOLUTION CF THE N79-23560 

EARLY EAETH: THE ROLE OF HAJCF IMEACTS 

(RASA) 36 p HC A03/HF A01 CSCL 08G 

Unclas 

G3/46 25920 

Crustal Evolution of the 
Early Earth: The Role of 
Major Impacts 


Herbert Frey 


MARCH 1979 


National Aeronautics and 
Space Administration 

Goddard Space Right Center 

Greenbelt, Maryland 20771 




CRUSTAL EVOLUTION OF THE EARLY EARTH: 
THE ROLE OF MAJOR IMPACTS 


Herbert Frey 

Geophysics Branch 
NASA Goddard Space Flight Center 
Greenbelt, Maryland 20771 

and 


Astronomy Program 
University of Maryland 
College Park, Maryland 20742 


Submitted to Precambrian Research 


ABSTRACT 


Impact cratering was an important— even dominant— process affecting the 
crustal evolution of the small terrestrial planets. The fundamental highlands/ 
maria dichotomy of the Moon's surface can be traced to a late heavy bombard- 
ment by basin-forming, asteroid-sized bodies which produced not only a 
topographic division in the lunar crust also localized the later eruptions 
of mare basalts. Major impact basins with diameters in excess of 200 km are 
recognized throughout the inner solar system from Mars to Mercury. Similar 
craters must have formed on the Earth prior to 4 billion years ago, and the 
minimum number of such basin-forming impacts can be calculated by scaling 
from the observed (minimum) number preserved on the Moon. When allowance 
is made for differences in impact velocity, gravitational cross-section and 
the effects of gravity on crater diameter, it is found that at least 50% of 
a presumed global sialic crust would have been converted into impact basins 
by 4 billion years ago. Among the effects resulting from the impact of an 
asteroidal object on the early crust were: (a) establishment of a topographic 
dichotomy of 3-4 km (after isostatic adjustment), (b) pressure-release 
partial melting of the upper mantle and rapid flooding of the basin floor 
by basalt, and (c) enhancement of thermal gradients in the sub-basin litho- 
sphere and upper asthenosphere. Comparative planetary data such as impact 
scaling can be used as important constraints on models of the early terres- 
trial crust. For example, the topography resulting from impact bombardment 
produced discrete oceans and dry land by 4 billion years ago, making unreason- 
able models of a globe-encircling ocean on the Earth after that time. 
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I NTRODUCT I ON 

The understanding of the present nature of the continental crust of the 
Earth requires knowledge of the evolution of the crust through time, as well 
as Information on the composition and structure of the original crust. While 
the recent evolution of the lithosphere may be described reasonably accurately, 
the accuracy of that description decreases rapidly the further into the past 
we "look. In particular, the scarcity of rocks on the Earth older than 
<3.5 billion years makes models of the earliest Archaean based only on 
terrestrial data highly speculative at best. Yet this is the very time when 
so many of the processes which today dominate the Earth's development were 
becoming established. Those processes have been sufficiently active that 
older rocks were obliterated as newer ones were produced. It is extremely 
unlikely that the clues necessary to understand the earliest Earth will be 
found only on this planet. 

But important constraints on modeiing the early Archaean are available. 

The smaller planets are generally less active and have experienced more 
sluggish evolution than the Earth; they have retained a greater fraction of 
the oldest rocks on their surfaces. In general evolutionary activity and 
complexity increases with size of the planet, because of the greater thermal 
reservoir available for those bodies. When allowance is made for differences 
in size, mass, bulk composition and location in the solar systems, studying 
the crustal evolution of the smaller terrestrial planets offers important 
information about the crustal evolution of the Earth as well (Frey and Lowman, 
1978). 
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This Is particularly true for early cratering effects. The Moon, Mercury, 
Mars and Venus all retain, to some degree, ancient cratered terrains. The 
tremendous geologic and meteorologic activity of the Earth erase traces of 
most craters within a few million years, and there is no evidence for any 
large craters dating from the early Archaean (Grieve and Robertson, 1978). 

Yet cratering was a fundamental, even dominant process throughout the inner 
solar system during the early stages of crustal formation. In this paper the 
role of major impact basins (such as those which formed on the Moon before 
4 billion years ago) is examined to determine the effects of such impacts 
on the early crustal evolution of the Earth. 

Specifically we address the following fundamental problem: what is the 
origin of the Earth's fundamental crustal dichotomy of low density continental 
and high density oceanic crust? Is this in any way related to the super- 
ficially similar highlands/maria crustal dichotomies of the Moon, Mercury 
and Mars? The lunar two-fold crustal division can be traced to late heavy 
bombardment by basin-forming, asteroid-sized objects. These impacts not only 
produced a topographic dichotomy, but localized the later eruptions of mare 
basalts which produced a compositltTJl dichotomy in the lunar crust. Is a 
similar series of events likely for the early Earth? 

The formation of the Earth's original crustal dichotomy is obscure. 

Today the crustal dichotomy is maintained by plate tectonic processes, but 
these same processes were not likely to have produced the original separation 
of high and low density crust. Comparative planetary data (Lowman, 1976) as 
well as terrestrial data (Shaw, 1972, 1976; Hargraves, 1976) suggest that 
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orlglnally the Earth was covered by a thin, low density, sialic crust. If 
true, plate tectonic processes could not have converted some 60% of that crust 
Into high density "oceanic" crust because of the difficulty of subducting con- 
tinental lithosphere once It has formed (McKenzie, 1973; Molnar & Gray, 1979). 
This especially true for the early Earth when the lithosphere should have 
been thinner. It can be argued (Frey, 1977a) that the crustal dichotomy must 
have been established first ; that Is, subductable crust must have been present 
before large scale horizontal displacement of the lithosphere was possible. 

But plate tectonic processes may have been active for much of the Earth's 
history (Burke et al., 1976; McKenzie and Weiss, 1975). Therefore It becomes 
desirable to look at the crustal evolution of the earliest Earth for the origin 
of the crustal dichotomy; that is, before 4 billion years ago, when cratering 
was still a very significant mechanism of crustal modification. 

BASIN-FORMING IMPACT CRATERING 

Large impact basins with diameters greater than 200 km are recognized 
throughout the Inner solar system, being well displayed In the ancient highland 
crusts of the Moon, Mercury and Mars (Wood and Head, 1976; Frey and Lowry, 1979). 
It would be impossible for the Earth to escape this bombardment which produced 
the lunar Imbrium and Oriental e basins before terminating about 4 billion years 
ago. This remains true regardless of whether the bombardment was part of the 
late accretlonal sweep-up (Wetherill, 1976) or part of a special flux of 
objects (Chapman, 1976). Dynamically plausible orbits exist for both cases 
(Wetherill, 1975, 1976). In this paper we do not consider whether the period 
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of basin-forming Impacts was a discrete "event" In the Inner solar system, 
but simply calculate the cumulative effects of such a bombardment until Its 
termination 4 billion years ago. Because these calculations are scaled from 
observed lunar Impact basins, they represent a minimum estimate of the total 
effects on the Earth. Such a scaling Is reasonable because the Earth and Moon 
would see approximately the same Incoming flux for most likely Impacting 
bodies from outside the Earth-Moon system. 

Consider objects In Apollo asteroid-type orbits, which are reasonable 
modern analogs of the bodies responsible for the late heavy bombardment of 
the Moon (Wetherlll, 1975). These objects approach the Earth-Moon system 
with a relative velocity between 15 and 20 km/sec (5p1k, 1966). Impact veloci- 
ties (Figure la) at the surface of the Earth will range from 18.7-22.9 km/sec, 
and at the lunar surface the objects will hit at 15.2 to 20.1 km/sec. The 
15-20% greater impact velocity at the surface of the Earth Is due to the 
greater gravitational acceleration of the Earth. Gravitational field differ- 
ences have additional effects, as described below. 

Faster impacts produce larger impacts on any planet. Impact velocity is 
generally related to crater diameter through an energy scaling law 


D * 


CE k = C (J*nV 2 ) 


k 


( 1 ) 


where D Is the crater diameter, E the energy released by the impact and taken 
to be equal to the kinetic energy of the Impacting body, C and k are constants. 
C is not well determined, but, if the same for the Earth and Moon, then the 
ratio of diameters of impact basins formed on the Earth to those formed on 


the Moon by objects of equal mass can be written 



APPROACH VELOCITY APPROACH VELOCITY 
(km/sec) (km/sec) 



Where and V, are the Impact velocities at the surface of the Earth and 
Moon, respectively. The constant k Is also uncertain. Values between 1/3 
and 1/4 have been suggested, but most authors adopt 1/3.3 or 1/3.4 (see 
Hartmann, 1965, 1977; Dence et al., 1977; Dence, 1978). In reality, the 
ratio above Is not very sensitive to k, as shown In Figure lb (Frey, 1977b). 
Here 1/3.3 Is used. For approach velocities of 15-20 km/sec, craters on the 
Earth might be expected to be some 11-15% larger than those formed on the 
Moon by Identical objects. In fact the basins formed on the Earth may be 
somewhat smaller than expected due to greater gravitational confinement of 
the cavity on larger planets. This can be represented by a gravitational 
acceleration ratio (e.g., Hartmann, 1977) 



The value of z Is uncertain. Hartmann (1977) adopted 0.2, but Dence et al. 
(1977) used a value of about 0.1. Below we treat three cases, retaining z 
as a variable yet to be determined. 

Therefore the diameter ratio of basins on the Earth and Moon Is 



Here we are concerned with the area of impact basins formed on the Earth com 
pared with the area of lunar impact basins. The ratio (D^/Dp) 2 Is plotted 
in the bottom portion of Figure 2a as a function of approach velocity, for 
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z « 0, 0.1 and 0.2. At 15 km/sec approach velocity the area ratio of basins 
ranges from 1.29 (z ■ 0.0) to 0.63 (z * 0.2). At 20 km/sec, basins on the 
Earth would be 1.17 (z ■ 0.0) to 0.57 (z ■ 0.2) times as large as those 
formed on the Moon by Identical objects. This clearly shows the Importance 
of the gravitational acceleration term z: for z * 0.2, terrestrial basins 
are much s maller than corresponding lunar basins despite the greater energy 
Involved In their formation. 

One other Important effect must be considered. Because of Its larger 
gravitational cross section, the Earth will collect many more of the Incoming 
bodies. This results In a much larger number of Impacts per unit area of all 
sizes. This In fact Is the dominant effect In these scaling relations, as 
shown elsewhere (Frey, 1977b). The gravitational radius 

i • 

Rg = 1 + _esc (5) 

where Rg is the gravitational radius, R the physical radius, V e$c the escape 
velocity, and the approach velocity. The ratio of the gravitational cross 
sections of the Earth and Moon Is shown In Figure 2a (times 10), as is the 
ratio normalized by the physical cross-section ratio of the two planets. 

This latter term Is the relevant one: 

( 6 ) 

as It shows, as a function of approach velocity, the ratio of the number of 
Impacts per unit area that will occur on each body. For approach velocities 
between 15-20 km/sec, the above term ranges from 1.52 to 1.30, decreasing with 
Increasing velocity. That Is, over the approach velocities of Interest, the 
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Earth collects 1.3 to 1.5 times as many craters per unit area as does the Moon. 

The effects of basin diameter and basin density may be combined into a 
single conversion factor 





, t» t 

(V r j > 


(7) 


This Is shown In F1gu,*e 2b, as a function of approach velocity, for three 
values of z. In this form It Is possible to calculate the equivalent area of 
basins formed on the Earth (In units of the surface area of the Earth) directly 
from a knowledge of the total area of lunar basins, normalized by the lunar 
surface area. This Is referred to as "equivalent basin area" and is denotv 
by a ©or ajj . Note that £ Is not the total surface area of a planet covets.: by 
basins, but the total area of all basins (normalized by surface area). In 
determining the total surface area covered by basins, allowance for overlap 
must be made, as described below. 

In principle all that is needed is knowledge of a^ , the equivalent area 
of all lunar basins In units of the Moon's surface area. Figure 3 is a histo- 
gram of all lunar basins larger than 200 km (Frey, 1977b). The total number 
of such basins is larger than that given by Wood and Head (1976), who restricted 
their count to multi ringed and peaked ringed basins. The histogram above does 
not Include Mare Gargantuan, and so should be considered a minimum for the 
number of lunar basins. The total area of these Is equivalent to 0.40 , 

where Is the surface area of the Moon. If Mare Gargantuan (Wood and Head, 
1976) is included the equivalent basin area becomes 0.50 A j . It is not unreason 
able that burled basins mcy raise this figure even higher, perhaps to 0.60 Aj> . 
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If the observed number of lunar basins are scaled to estimate the number 

of basins forming on the Earth, the terrestrial basins will be an under-estimate 

of the number that probably formed. The reasons for this are two. First, as 

discussed above, the actual number of basins that formed on the Moon Is almost 

certainly larger than the number preserved. Secondly, there Is the distinct 

possibility that the Earth sustained even larger impacts than those seen on the 

Moon. As described elsewhere (Frey and Lowry, 1979) the lunar basins, when 

plotted on a log (cumulative number/ km )-log (diameter) graph, closely follow 

a line with slope = -2 (for D£ 300 km). If we assume a D distribution of 

basins on the Earth as well, then it is possible for a number of very large 

basins to have formed before 4 billion years ago. Table I shows this. If 

there are 3 basins on the Moon with diameter equal to or greater than 1000 km, 

then there should have been some 56 basins formed on the Earth larger than or 

equal to 1000 km diameter. This in turn "predicts" 6 basins larger than 3000 

_o 

km, 2 larger than 5000 km, and 1 larger than 7000 km, according to a D law. 

If there are 4 lunar basins larger than 1000 km, then the Earth might have had 
8 larger than 3000 km and 3 larger than 5000 km. These are approximate values, 
and of course depend on the value of C® adopted (that is, on the ratio of 
basin areas on Earth and Moon and the number per unit area ratio for the Earth 

and Moon). Figure 4 shows those basins which might be ignored unless this 

_2 

effect is taken into account. The solid line is a D fit through D=300 km 
on the log (cumulative number) vs log (diameter) plot. The long dashed line 

_p 

is the corresponding D distribution expected tor the Earth, for = 1.4 
(see below). This corresponds to a value of z = 0.1 at V w = 15 km/sec. The 
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lunar distribution cuts off at 2100 km; there could be 6 basins larger than this 
formed on the Earth as shown by the extrapolation of the Earth D" 2 line (short 
dashes). The total area of these "extra" basins Is some 1.23 x 10 8 km 2 , or 
roughly 0.24 A$ , where A^ is the surface area of the Earth. This compares 
with the total area of about 0.55 A^ for all the basins with Dl 2100 km that 
are expected to have formed on the Earth based on the number seen on the Moon. 

It Is clear that these "extra" basins represent a significant contribution to 
the total equivalent area of terrestrial basins. 

The above effects are combined in Figure 5. Three cases are shown, for 
three values of the approach velocity (15.0, 17.5 and 20.5 km/sec). Approach 
velocity contributes both to impact velocity (and therefore the individual 
basin area ratio [.D^/Dj^ 2 ) and to the gravitational cross-section ratio 
(Rg® /Rg^ ) 2 . Consider the middle panel, for * 17.5 km/sec. In this (and 
the two adjacent panels) the horizontal axis is a^ , the equivalent lunar 
basin area (total area of all lunar basins divided by surface area of the Moon) 
As described above, this value is at least 0.4 the lunar surface area (although 
basin overlap reduces the actual surface area covered by basins to about 0.32 
Aj ). The vertical axis is the equivalent basin area on the Earth a$ . That 
is, a^ is the total area of all basins produced on the Earth by scaling from 
lunar basins. The minimum value of this is given by 

a © = a j> 

where C$ is the conversion factor shown in Figure 2b. C$ of course depends 
on the gravitational factor z. For z = 0, the diameters of terrestrial basins 
are much larger than those formed on the Moon by identical objects (see Figure 
2a). Thus is a maximum of this value (1.69). For Z = 0.2, is reduced 
to 0.83, because terrestrial basins have only 60 % as much area as the corres- 
ponding lunar basins (see Figure 2a). 


15.0 km/sec [ V = 17.5 km/sec | V = 20.0 km/sec 
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Equation 8 above assumes no terrestrial basins larger than the largest 
lunar basin formed on the Earth. For this case, the calculated value of a^ as 
a function of a^ is given by the solid line in Figure 5. However, "extra" 

-2 

basins may have formed on the Earth due to the large diameter fact in the D 

distribution (Figure 4). The elevated dashed line in Figure 5 accounts for 

these extra basins, whose number and size depend on the approach velocity and 

the assumed value of z. The true a^ probably lies somewhere between the solid 

and dashed lines, for any given value of aj> . That is, a^ is probably larger 

than that given by the lunar size distribution alone, but is probably less than 

-2 

the maximum possible using a pure D law for terrestrial basins. 

In fact, the value of a^ is not enough. This only represents the total 
area of basins formed on the Earth. A more significant quantity is the total 
fraction of the Earth's surface converted into basins, A. From the probability 
of each new impact basin forming in an uncratered area, it can be shown (Frey, 
1977b) that 

A = 1 -e' a » (9) 

The scale for A, In units of the Earth's surface area, Is shown on the right. 

The scale is compressed at the bottom as would be expected: as the total area 
(or number) of basins Increases, it becomes harder and harder to find uncratered 
areas, so the total fraction of the Earth's surface not cratered decreases more 
slowly. That is, most impacts occur in areas already cratered. 

For a middle velocity of 17.5 km/sec and assuming z * 0.1, = 1.18. 

That is, the area of basins on the Earth is at least 18% greater than those 
on the Moon, both normalized by the planetary surface area. Again, the probable 
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existence of "extra" basins suggests that the effective C$ Is greater than 
this. If there are lunar basins corresponding to 40* the lunar surface area, 
the equivalent area of basins on the Earth Is 0.47-0.64 Ag , which corresponds 
to 38*-47* of the Earth's surface converted Into basins when allowance Is made 
for overlap. A more reasonable value for a^ might be 0.5Aj , to allow for 
uncounted lunar basins (such as Mare Gargantuan). Then 1.18 (z > 0.1 at 
17.5 km/sec) Implies * 0.59-0.76 A^ and A = 0.45-0.53A$ . 

In general a lower approach velocity Increases the total basin area on 
the Earth as a function of a^ and z. This Is because lower velocity of 
approach Increases the ratio of Impact velocities on the Earth to those on the 
Moon (see Figure 1), making the diameters (for any given z) generally larger 
on the Earth (Figure 2a). Lower approach velocities also Increase the ratio of 
gravitational cross-sections (Earth/Moon), meaning the Earth sustains more 
impacts/unit area because Its greater gravity can pull In more slow moving 
objects. Therefore Increases sharply with decreasing approach velocity 
(Figure 2b). Consider s 15 km/sec, z = 0.1 again. Then for a^ s 0.4 A j , 
a^ = 0.54-0.71 A$ (allowing for extra basins) and A-0.42-0.51 A 0 . For a^ * 
0.5 Aj , these Increase to a$ = 0.68-0.85 A$ and A = 0.49-057 A^ . 

Therefore It would seem that the observed number of lunar basins plus 
those "extra" basins suggested by the D law yield a basin equivalent area of 
at least 0.70 A^ , which In turn means at least 50* of the Earth's surface was 
converted Into basins by 4 billion years ago. This value of "at least 50*" 
seems warranted considering that the number of lunar basins observed Is truly 
a minimum, and agrees with a previous estimate made using simpler arguments 
(Frey, 1977b). Note the highest value suggested by Figure 4 (V^ = 15 km/sec, 
z s 0 and a^ = 0.6) has about 70* of the Earth surface covered by basins. 
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EFFECTS OF BASIN-FORMING IMPACTS 

There are three major effects that mult from the formation of a large 
Impact basin: (a) creation of a large but shallow crater, (b) pressure-release 
partial melting at depth and subsequent flooding of the basin by basaltic lavas, 
and (c) alteration of the geothermal gradients below the basin (Frey, 1977a). 

All of these effects depend on the depth of the crater formed. It Is not 
obvious that either lunar, mercurlan or terrestrial depth (d)-dlameter (D) 
relations apply to the case of very large Impact basins (diameter > 1000 km). 

All three planets show a distribution of depths that can be described by two 
straight lines. For the Moon (Pike, 1974) 

d * 0.196D 1,010 D 10km (10a) 

d = 1.044D°' 301 D 10km (10b) 

Mercury craters follow this broken distribution to within 10% for fresh craters 
(Mai In and Dzurzln, 1978). The transition to a lower slope probably Involves 
modification of a the crater form during Impact. Terrestrial craters change 
from simple to complex (for those formed In crystalline rock) at about 3-4 km 
diameter (Grieve and Robertson, 1978) while lunar and mercurlan craters change 
slope at about 10 km (Mai In and Dzurzln, 1978). Grieve and Robertson (1978) 
give for the Earth 

d * 0.3260®*^ (simply i n crystalline rock) (11a) 

d * 0.52D 0,1 ® 9 (complex, In crystalline rock) (lib) 

but point out the large scatter in the d-D relation for complex craters makes 
the slope so uncertain that It could easily be^O^, as is the case for lunar 
craters. 
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The lunar curve for large craters, If extrapolated Into the diameter range 
for large Impact basins, suggests a 1000 Ion basin would be^8 km deep, the 
terrestrial curve gives a value of ~2 km deep. These very shallow depths are 
undoubtably due to crater modification by a variety of processes Including 
Isostatic adjustment, wall failure and si imping both during and following 
crater formation. Extrapolation of the simple crater curves, which might be 
thought of as representing a transient cavity or "Instantaneous crater form," 
yield incredible depths of 210 km (from lunar curve) and 74 km (from terres- 
trial curve). These depths transcend the crust and, In the extreme, the litho- 
sphere. Under these circumstances It Is very unlikely that normal cratering 
forms and processes apply. The equation of state of the target material 
changes over the presumed depth to which a shock wave might penetrate; propa- 
gation of a shock wave will no longer be spherical and the resulting cavity 
may be significantly more shallow that for craters formed In layers which are 
thick and homogenous compared to crater dimensions. 

Here we are dealing with only the first order effects of major Impacts; 
details of crater forms are beyond the scope of this paper. Below it is 
assumed that initial crater depths were shallow: a 1000 km wide basin is 
assumed to have been M3 ion deep at the outset (Baldwin, 1963). All of the 
effects described below are significantly enhanced If the basin was Initially 
much deeper. The descriptions Included here should therefore be viewed as a 
conservative "worst case." 
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To describe even the mechanical effects of a major Impact, a model for the 
early Earth lithosphere Is required. This becomes even more crucial for dis- 
cussion of thermal effects. For reasons discussed earlier It Is assumed that 
the Earth originally had a global crust of Intermediate composition (Lowman, 

1976; Shaw, 1972). This crust Is represented as andesite (Lowman, 1978) and 
assumed to have been 20 km thick. The crust overlies a pyrollte mantle below 
which Is perldotlte (Frey, 1977a, c). Consider a late-arriving asteroldal 
object which Impacted this crust roughly 4 billion years ago and formed a 
basin 1000 km wide. 

Initially the crater was ~13 km deep, but, In the thin lithosphere of the 
Earth this basin would soon isostatlcally adjust upward. The basin floor 
would rise over several thousand years, and leave a relief between basin rim 
and floor of 3-4 km (Frey, 1977a). This value is relatively Insensitive to 
the Initial depth of the crater, but is quite dependent on the assumed thickness 
of the low density crust. Decreasing the crustal thickness decreases the total 
adjusted relief, as a larger percentage of the excavated material is then of 
higher density. 

Formation of such a crater produced a pressure drop at the surface. 

Removal of the overburden means that pressure- temperature relations at depth 
were altered In the direction that favors melting at more shallow depths. 

Once solid material near Its pressure-melting temperature would have partially 
melted when the pressure was reduced. For a pyrolite mantle the partial melt 
would have been basalt; the liquid would have been less dense than the surrounding 
solid but highly fractured rock and should have risen rapidly through the 
thinned and broken lithosphere (Frey, 1977b). Basin filling by basalts would 
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have begun almost Immediately following Impact. Note the above has not taken 
account of any extra heat deposited at depth due to the Impact Itself (Safronov , 
1978). The appearance of basaltic material In the topographic lows of major 
Impact basins produced a compositional dichotomy on the early Earth akin to 
that of the Moon. By the time the bombardment had come to an end *4 billion 
years ago, at least 50% of an original global sialic (continental) crust had 
been converted Into terrestrial maria. These maria were the original oceanic 
basins of the Earth, as described elsewhere (Frey, 1977b). Thus the funda- 
mental crustal dichotomy of the Earth was established very early In Earth 
history, making possible early plate-tectonic processes requiring subductable 
crust. 

Note that If the Impact was initially deeper than described above, the 
flooding of the basin would have been even more rapid. Deeper excavation pro- 
duces a greater pressure drop leading te iiare complete melting at shallower 
depths. The liquid then had even less distance to rise to erupt onto the basin 
floor. In the extreme limit of depths^ 7C km the impact would have penetrated 
Into already molten regions and basaltic lava would have been immediately 
exposed. 

It has been suggested that some highly altered remnant of this original 
mare-type volcanlsm may still exist on the Earth as the greenstone belts found 
in Archaean terrain (Green, 1972; Glickson, 1976). Many ultrabasic Archaean 
rocks seemed to have formed In the presence of unusually steep thermal melting 
gradients (Green, 1975), such as might be associated with large impacts (see 
below) and the corresponding pressure change such impacts produced. No known 
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rocks on the Earth data from 4 billion years ago, but greenstones with ages of 

3.5 billion years old do exist. At best such rocks would be late eruptions. 

Other possible modes of origin are also likely (see Wc’blen and Schulz, 1978; 

Burke et al., 1976; Gllckson, 1976; Green, 1975). 

The question of early geothermal gradients and their variation due to 

Impact Is an Important problem. In order to examine the first order effects a 

simple thermal history was calculated as follows. The Earth was assumed 

Initially molten, but differentiated Into a 20 km thick andesitic (or dlorlte) 

-3 

crust (density 2.7 gm cm ), an upper mantle of pyrollte (density 3.3 from 20 
to 300 km depth) and a lower mantle of perldotlte (density 3.5 from 300 to 500 
km depth). Specification of these rock types determines the radioactive heat 
contributions from U 22 ®, U 252 , Th 222 and K^®. The one dimensional equation for 
thermal conduction Is 



where K Is the mean thermal conductivity (here taken to be 0.008 cal cm’*sec”* °C”*), 
Cp Is the specific heat at constant pressure (0.25 cal gm"* °C"*), p Is the 
density, T the temperature, t Is time and H(Z,t) Is the rate of heat generation 
per unit volume due to radioactivity. For this simple case consider only the 
outer 500 km of the Earth, divided Into layers 10 km thick. Assume the lower 
boundary has fixed temperature. The surface temperature Is set equal to zero 
and the equation solved by finite differences In time steps of 1 million years 
(a requirement fixed by the stability crlterlal for the finite difference tech- 
nique). After 500 million years (l.e. roughly 4.1 billion years ago) the 
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evolution Is interrupted to allow for a 1000 km wide, 13 km deep Impact. The 
newly exposed surface layer temperature is set equal to zero, the pressure and 
melting temperatures at depth recalculated to allow for the excavation, and the 
evolution allowed to continue for another 500 million years. Details are 
available elsewhere (Frey, 1977c). 

Some of the results are shown in Figure 6. Here depth to various Isotherms 
is plotted as a function of time since origin. The line labeled D refers to 
the depth at which the temperature exceeds the melting temperature at that 
depth; this is taken to be the lithosphere/asthenosphere boundary. Solid lines 
show the effect of the impact (labeled B, at 500 million years after origin) 
on the Isotherms. Dashed lines show what would be the case If no Impact had 
occurred; this may be taken as the thermal profile for adjacent highlands. 

Initially the isotherms deepen with time as the molten model cools, but 
then begin to rise slightly after 150 million years as radioactive heating 
becomes important. The lithosphere deepens to about 75 km after 100 million 
years, but then thins to about 64 km just prior to the basin-forming impact. 

Had no impact occurred, the isotherms would remain more or less stable in this 
simple model, especially over the outer few hundred kilometers. 

The effect of the impact is obvious. All isotherms except the melting 
line (0) drop rapidly. This is a consequence of penetration to deeper, hotter 
layers which cool rapidly by radiation, causing the loss of heat at depths to 
Increase. The depth to melting (strictly speaking, partial melting) Initially 
rises several km immediately after impact. This is due to the pressure- 
released melting associated with excavation of the basin. Melting occurs 
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closer to the surface than before. The lithosphere Is therefore vhlnner, 
especially due to the Impact having stripped off at least 13 km of crust* as 
shown In Figure 6. Following the “D" Isotherm also shows how the thickness 
of the lithosphere changes with time. All the Isotherms Increase with depth 
as the model cools; 500 ml 11 Ion years after the Impact "D" has deepened to 79 km. 
But the excavation of the crater removed 13 km of crust so the true thickness 
of the lithosphere Is only 66 km, roughly equal to Its pre-impact value. That 

O 

Is, In this simple conductivity model It takes some 10 years for the litho- 
sphere below the basin to cool enough to thicken to the highland lithospheric 
thickness. 

This time to cool will obviously depend on the depth of the impact; 
deeper penetration to still hotter layers results in more rapid cooling and 
thickening. For a 10 km deep impact, the time to rethicken to highland litho- 
spheric values is 100 million years; if the impact was 50 km deep then 220 
million years are required (Frey, 1977c). Qualitatively these results are in 
agreement with a similar case done for the iioon (Arkani-Hamed, 1974), although 
subbasin lithospheric thickening is some 2.5 times slower than In the lunar 
case. The major point here is that the subbasin lithosphere remains thin 
(mostly due to impact excavation) compared with the highland lithosphere for 

Q 

some 10 years following impact. 

Formation of such a large basin would also have altered geothermal gradi- 
ents. Exposure of hot deep layers to rapid cooling by radiation created very 
high gradients across this new surface. Gradients at depth were also affected, 
as shown in Figure 7 (Frey, 1977a, c). Here gradients averaged over 20 km thick 
layers a e shown as a function of time before and after the basin-forming impact 
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(solid lines). The no-impact case is shown by dashed lines. At all depths 
thermal gradients are enhanced. In the 30-50 km layer (13-37 km below the 
new basin surface after impact) the change is from 18°C/km before to 23°C/km 
immediately after Impact. At greater depths the effect is somewhat less, but 
it is clear that thermal gradients were enhanced by about 20% below the litho- 
sphere by the impact. Note also that this steepening of the gradients persists 
over time. Had the impact been deeper than the 13 km used here the enhance- 
ment would have been even greater due to exposure of even hotter layers to 
rapid cooling. 

The above simple thermal history has not included possible convection on 
the Earth, and is therefore unrealistic. McKenzie and Weiss (1975) have 
shown that some form of mantle convection has most likely always been part of 
the Earth's solution. Although not included explicitly, the qualitative 
effects of the basin-forming impact on early convection are obvious. Steepening 
the thermal gradients by 20% stirred up convective transport below the basin 
by 20%. This more vigorous convection operated below a thinned (by at least 
25%) and badly fractured lithosphere. These were conditions favorable for 
rapid rifting or break up of the subbasin-lithosphere. It is not unreasonable 
that early plate formation and plate motion may have rapidly followed the 
formation of major impact basins and their flooding by basalt. That flooding 
provided subductable crust at the surface. The badly fractured lithosphere 
could easily have been broken into a large number of small microplates 
(Frey, 1977a, c) which might be expected to have moved more rapidly than the 
large plates of today. McKenzie and Weiss ( 1975 ) and Burke et al. (1976) 
have also argued for small plates and/or rapid motions in the Archaean. 


- 21 - 


The above discussion does not Include the possibility of significant 
amounts of heat being deposited within the lithosphere by the Impact Itself 
(Safronov, 1978). This Is a difficult problem to assess quantitatively, 
because for Impacts of the size described here there are enormous uncertainties 
ir. energy scaling, energy partition and the depth dependence of the thermal 
energy deposited. However, It seems that most high temperature effects due 
to Impacts are concentrated near the surface (Grieve et al., 1977). It may 
be that the near surface thermal gradients were overturned, but this should 
have been a short-term effect compared with others discussed here. There is 
not likely to have been a significant effect at tne depths of the lithosphere/ 
asthenosphere boundary unless the impact penetrated this far; in this case 
the incremental heat added at these levels (60-70 km depth) would not 
significantly add to the thermal energy of the layer (see Figure 6). 

DISCUSSION 

It is clear that the Earth, like the other terrestrial planets, must have 
experienced severe Impact bombardment early in its m story by basin-forming, 
asteroid-sized objects. A conservative scaling from the observed lunar basins 
indicates that the terrestrial bombardnent was enough to convert at least 
half of a global sialic crust into basins some 4 km deep. For any reasonable 
thermal profile for the early Earth the lithosphere would have been thinner 
than today. Eruption of basaltic lavas onto the basin floor is expected due to 
pressure-release melting below the impact. A crustal dichotomy like the lunar 
highlands/maria division and superficially like the modern continents/ocean 
basins was established early in the history of the Earth, by 4 billion years 


ago. 
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This is the general picture provided by comparative planetary evolution. 

The details require more study. But the above-described scenario is more 
realistic than many stories written about the early Earth because it is based 
on observational data. Too often models of the Earth are constructed based on 
terrestrial data alone. These are generally unsatisfactory for the earliest 
evolution of the Earth because of the lack of terrestrial data from the early 
Archaean. Many such models have completely ignored the effects of impacts on 
the early crust. One example is discussed below. 

Hargraves (1976) has suggested a model for the primitive Earth in which 
a global sea of some 2 km depth overlies a primordial sialic crust. In his 
model continental crust thickens with time, finally emerging above sea level 
quite late in Earth history at 1.7 to 2.3 billion years ago, or, in his pre- 
ferred version, not until the late Precambrian. Comparative planetary data 
would suggest that such an emergence of dry continental crust occurred much 
earlier, by no later than 4 billion years ago, and that it was due to the same 
impact bombardment described above (Frey, 1979). 

If the early Earth was originally covered by a 2 km deep global sea, the 
intense cratering by basin-forming objects would have caused sea level to pro- 
gressively lower. This is shown in Figure 8. If we assume an average adjusted 
basin depth of ^4 km, then the volume of basins increases linearly (as shown 
by the solid line) as the percentage of the Earth's surface converted into 
basins increases. The volume of water on the surface of the Earth minus 
the cumulative volume of basins formed decreases at the same rate, as shown by 
the large dashes. The two curves cross at 0.25 of the Earth's surface converted 
to basins; at this point the cumulative volume of basins equals half the original 
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volume of water and sea level has dropped from 2 to 1 km average depth (short 
dashed line). Impacts which produced basins 1000 km across also produced rim 
relief on those basins In excess of 1 km, if lunar cratering Is any guide 
(Pike, 1967). Therefore the highest basin rims were above sea level when 
basins covered only 25% of the Earth's surface. As an Increasing percentage 
of the surface was converted Into 4 km deep basins, water continued to drain 
Into these sinks, and sea level continued to drop. Figure 9 is a schematic 
plot that represents this progressive emergence of dry "continental" or high- 
land crust. By the time the bombardment came to an end 4 billion years ago, 
at least 50% of the Earth's surface had been converted to basin topography; 
sea level had lowered to about the original crustal radius of the Earth leaving 
nearly half of the surface (at least slightly) above sea level. 

The above says nothing about the degassing history of the Earth. It is 
not clear that a global ocean did exist on the primordial Earth. If it did, 
the topography produced by the basin-forming impact bombardment prior to 4 
billion years ago reduced sea level to a point where a large percentage of the 
Earth's "continental" crust stood dry by the end of the bombardment. If 
degassing was less catastrophic, the already existing basins would have provided 
natural sinks into which late accumulating water would drain. It seems clear, 
therefore, that dry continental crust has been present on the Earth from at 
least 4 billion years ago, and did not remain submerged for much longer as 
Hargraves (1976) suggested. 

The implications of an early "emergence" of continents have been discussed 
elsewhere (Prey, 1979). The point to be made here is the value of comparative 
planetary data in constraining models of the early Archaean. Scenarios which 
ignore the obviously important effects of early cratering are not likely to be 
realistic. The late heavy bombardment was a dominant process affecting the 
crustal evolution of the smal 1 terrestrial planets; it was also significant 
in the crustal development of the Earth. 


GLOBAL OCEAN INITIALLY 2km DEEP 
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CONCLUSIONS 

The lack of terrestrial data from 3-4 billion years ago severely restricts 
the accuracy of models of the early Archaean based on this information alone. 

The more sluggish crustal evolution of the smaller terrestrial planets has 
preserved a more complete record of the earliest processes affecting these 
bodies. When allowances for differences in size, mass, bulk composition and 
location in the solar system are made, the data from the smaller planets can be 
used to constrain models of the early Earth. This is especially true for the 
dominant external process which affected the crusts of all terrestrial planets, 
the late heavy bombardment by asteroid-sized objects. These impacts occurred 
throughout the inner solar system, and terminated on the Moon about 4 billion 
years ago. The Earth could not have escaped such impacts, and they must have 
significantly modified the early crust of this planet. That crust was probably 
thin, intermediate In composition and global in extent, again based on compara- 
tive planetary analogy. 

An estimate of the effects of a late heavy bombardment on the early terres- 
trial crust can be made by scaling from the observed number of lunar Impact 
basins with diameters in excess of 200 km. This number is certainly a minimum 
and therefore any effects should be considered lower limits for the Earth. 
Craters forming on the Earth will be larger than those forming on the Moon due 
to identical objects unless gravity severely restricts the growth of the cavity. 
The ratio of impact basin area on the Earth to that on the Moon increases as 
approach velocity decreases , because the ratio of impact velocities increases 
as approach velocity decreases. A much more significant effect is that of the 
relative gravitational cross-sections of the Earth and Moon: for reasonable 
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approach velocities the Earth should have collected 1.3-1. 5 times as many Impacts 
per unit area as did the Moon. By the time the lunar bombardment came to an end 
roughly 4 billion years ago at least 50% of the Earth's original global sialic 
crust would have been converted Into large basins. Some 60 of these would have 
□ten more than 1000 km across. Furthermore, it Is possible that the Earth sus- 
tained a few very large impacts, which would have produced basins larger than 
those seen on the Moon. 

There were several immediate effects of these large-basin-forming Impacts. 
After isostatic adjustment typical basin floors would have been some 4 km below 
the cratered highland surface. Thus the bombardment created a topographic divi- 
sion with half of the Earth at higher elevations. These low lying basins would 
have been rapidly flooded by basalts produced at depth below the basin by 
pressure-release partial melting. Excavation of the crater provided the reduced 
pressure. Penetration to deeper, hotter layers caused these levels to cool 
rapidly and steepened the thermal gradients below the basin. This would have 
stirred up whatever convection existed in the asthenosphere. Combined with the 
highly fractured lithosphere, it is plausible that rapid rifting and breakup of 
that lithosphere produced a number of small and rapidly moving "microplates." 

If the degassing of the Earth was complete enough to have produced an early 
global ocean overlying the original crust, the topography which resulted from 
Impact bombarctnent lowered sea level such that dry continental masses existed 
before the bombardment was complete. The large impact basins provided natural 
sinks into which more and more water drained as an increasing percentage of the 
Earth's surface as converted into these basins. Models which suggest a global 
ocean persisting until the late Precambrian are unlikely. 
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The late heavy bombardment by basin-forming asteroldal-slzed objects which 
terminated roughly 4 billion years ago on the Moon was a dominant process affecting 
the crustal evolution of the smaller terrestrial planets; It was also significant 
In the crustal development of the Earth. That bombardment was responsible for 
the original ocean basins on this planet and may have triggered plate tectonic 
processes. These effects should not be Ignored In models of the primitive Earth, 
and can best be understood In terms of comparative planetary studies of the 
evolution of the terrestrial planets. 
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FIGURE CAPTIONS 

Figure 1: Impact velocity and relative crater diameter as a function of approach 

velocity, (a) Impact velocity of the surface of the Earth (©) and 
Moon & } and the ratio of these. Note that the ratio increases with 
decreasing approach velocity due to the greater effect of gravitational 
acceleration on slow-moving objects, (b) Basin diameter ratio versus 
approach velocity. Crater diameter is related to impact velocity 
through an energy scaling law. The ratio of diameters is not very 
sensitive to the exponent k, as shown by the two plotted curves. 

Figure 2: (a) Ratio of areas of basins formed on Earth to those formed on the 
Moon by identical objects for different values of a gravitational 
acceleration parameter (bottom). Larger values of Z correspond to 
greater confinement of crater cavity, resulting in terrestrial basins 
being smaller than lunar basins despite greater energy released at 
impact on the Earth, (top) Ratio of gravitational cross-section for 
the Earth and Moon (xlO) versus approach velocity, and ratio of gravi- 
tational cross-sections normalized by physical cross-section ratio. 

At lower approach velocities the Earth collects many more impacts per 
unit area than does the Moon, (b) Conversion factor versus 
approach velocity for different values of the gravitational accelera- 
tion parameter. This factor combines the ratio of basin areas and the 
gravitational cross-section ratio normalized by physical cross-section 
ratio (a) into a single term which allows calculation of the equiva- 
lent surface area of basins formed on the Earth from the total area of 
basins formed on the Moon. This yields a minimum total area of terres- 
trial basins (see text). 
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Figure 3: Histogram of lunar basins with diameter larger than or equal to 

200 km from Frey (1977b). All craters are counted without regard 
to specific morphologic subtypes. 

Figure 4: log (cumulative number)-log (diameter) plot for observed lunar 

basins (Figure 3) and calculated terrestrial basins. Solid line is 
_2 

D slope fitted through D=300 km. Long-dashed lines represent 
scaled terrestrial basins based on observed lunar basins with 
Di 2100 km. Short dashed line is extrapolation to possible addi- 
tional large diameter basins which might have formed on the Earth but 
did not on the Moon. 

Figure 5: Calculated area of basins formed on the Earth as a function of total 

lunar basin area. Both are normalized by planetary surface area 
(i.e., a^ = total area of lunar basins divided by lunar surface 
area). Each panel corresponds to a given velocity of approach (shown 
at the top). Three values of Z (0.0, 0.1 and 0.2) are shown. The 
solid line represents a^ calculated from (Figure 2b), that is, 
from extrapolation of only observed lunar basins. See text for 
details. Scale on right shows actual fraction of the Earth's surface 
covered by the total basin area given on the left scale. 

Figure 6: Depth to selected isotherms versus time for a simple, purely conduc- 

tive thermal history. "D" represents depth to (partial) melting, 
here taken to be the depth to base of the lithosphere. A 1000 km wide 
impact basin is shown as having occurred at 500 million years after 
origin. 
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Figure 7: Change in thermal gradient in 20 km thick layers versus time for the 
same thermal history model shown in Figure 6. After impact the 30- 
50 km deep layer is only 17-37 km below the new surface layer. All 
thermal gradients are enhanced by the impact which is introduced 

500 million years after orign; those nearest the impact surface are 

steepened the most. 

Figure 8: Change in sea level with increasing fraction of the Earth's surface 
covered by basin. After adjustment basins are assumed to be 4 km 
deep. Total volume of basins increases with time, causing sea lavel 
to decrease as water drains into the depressions. Global sea is 
assumed to be 2 km deep at the outset. 

Figure 9: Schematic diagram of progressive changes in sea level with increasing 
basin area. When 25% of the Earth's surface has been converted into 

basins 4 km deep, sea level has dropped to 1 km and the highest 

basin rims have already emerged from the sea. When 50% of the Earth 
is covered by basins, roughly half of the Earth's crust stands above 
sea level. 



